***US level regressions
set more off
use G:\paris\hurricane\state\growth_gdp_q_us, clear
sort year quarter
merge year quarter using G:\paris\hurricane\state\growth_pi_q_us
drop _merge
sort year quarter
merge year quarter using G:\paris\hurricane\revisions\hurr_us_q
*replace hurr=0 if hurr==.
replace hurr_c=0 if hurr_c==.
drop _merge
quietly tab year, gen(td)
quietly tab quarter, gen(qd)
sort year 
replace hurr=hurr/1000000
replace hurr_c=hurr_c/1000000
local i=1
while `i'<4{
	gen hurr_c`i'=hurr_c[_n-`i']
	local i=`i'+1
	}
*keep if year>=1948
*drop if year==1948 & quarter==1
replace growth_pi_q_us=growth_pi_q_us/4
replace growth_gdp_q_us=growth_gdp_q_us/4

drop if year>2005
capture log close
log using g:\paris\hurricane\revisions\national, append

reg   growth_pi_q_us initial_pi_q_us  hurr_c* td* qd* if year>=1970, vce(hc3)
outreg initial_pi_q_us hurr_c*    using g:\paris\hurricane\revisions\table6.txt, replace nolabel se coefastr nonotes bdec(3) nor 
reg   growth_pi_q_us initial_pi_q_us  hurr_c* td* qd* if year>=1948, vce(hc3)
outreg initial_pi_q_us hurr_c*    using g:\paris\hurricane\revisions\table6.txt, append nolabel se coefastr nonotes bdec(3) nor 
reg   growth_gdp_q_us initial_gdp_q_us  hurr_c* td* qd* if year>=1970, vce(hc3)
outreg initial_gdp_q_us hurr_c*    using g:\paris\hurricane\revisions\table6.txt, append nolabel se coefastr nonotes bdec(3) nor 
reg   growth_gdp_q_us initial_gdp_q_us  hurr_c* td* qd* if year>=1948, vce(hc3)
outreg initial_gdp_q_us hurr_c*    using g:\paris\hurricane\revisions\table6.txt, append nolabel se coefastr nonotes bdec(3) nor 

log close